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Abstract 

We compute QCD corrections to the processes e + e~ — > tihP and pp — ► iih by treating 
the Higgs boson as a parton which radiates off a heavy quark. This approximation is 
valid in the limits « Mf and Mh,Mt « y/s. The corrections increase the rate for 
pp — ► at the LHC by a factor of 1.2 to 1.5 over the lowest order rate for an intermediate 
mass Higgs boson, Mh < 180 GeV. The QCD corrections are small for e + e~ — ► tihP at 
y/s= 1 TeV. 



1 Introduction 



The search for the Higgs boson is one of the most important objectives of present and future 
colliders. A Higgs boson or some object like it is needed in order to give the W ± and Z 
gauge bosons their observed masses and to cancel the divergences which arise when radiative 
corrections to electroweak observables are computed. However, we have few clues as to the 
expected mass of the Higgs boson, which a priori is a free parameter of the theory. Direct 
experimental searches for the Standard Model Higgs boson at LEP and LEPII yield the limit, [|TJ 

M h > 77.5 GeV , (1) 

with no room for a very light Higgs boson. LEPII will eventually extend this limit to around 
Mh > 95 GeV. Above the LEPII limit and below the 2Z threshold is termed the intermediate 
mass region and is the most difficult Higgs mass region to probe experimentally. 

In the intermediate mass region, the Higgs boson decays predominantly to bb pairs. Since 
there is an overwhelming QCD background, it appears hopeless to discover the Higgs boson 
through the bb channel alone. The associated production of the Higgs boson could offer a tag 
through the leptonic decay of the associated W boson or top quark, 



pp -> Wh° 

pp -> tih° . (2) 

Both of these production mechanisms produce a relatively small number of events, making it 
important to have reliable predictions for the rates. 

The QCD radiative corrections to the process pp — > Wh° have been computed and increase 
the rate significantly, (~ 30 — 40% at the LHC) ||f§]. The QCD radiative corrections to the 
processes e + e~ — > tih° and pp — > tth° do not yet exist and are the subject of this paper. 

There is considerable expertise available concerning QCD corrections to Higgs production. 
The radiative corrections of O(a^) to gg — > h° involve two-loop diagrams and have been cal- 
culated with no approximations ||. They have also been computed in the limit in which 
Mh/Mt — > 0, where an effective Lagrangian can be used and the problem reduces to a one- loop 
calculation [0. In the intermediate mass region, one expects this to be a reasonable approxi- 
mation. The QCD corrections to the process gg — > h° can be conveniently described in terms 
of a K factor, 

K(u\ ~ a{PP - h ° X)NL ° W 

where \l is an arbitrary renormalization scale, which we take to be /i = Mh- Fig. [I] shows the K 
factor at the LHC computed exactly and in the Mh/Mt — > limit. The important point is that 
the Mh/Mt — > limit is extremely accurate and reproduces the exact result for the K factor to 
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Figure 1: Ratio of the radiatively corrected cross section to the lowest order result for pp — > h° 
at the LHC, y/s = 14 TeV, with \i = M h . This figure uses the 2— loop value of a a ((j) for both 
the lowest order and the NLO predictions. 

within 10% for all Mh < 1 TeV, [[| although the small Mh/M t limit does not accurately give 
the total cross section for M/, > M t . The reason the K factor for this process is so accurately 
computed in the small Mh/Mt limit is that a significant portion of the K factor comes from a 
constant rescaling of the lowest order result, along with a universal contribution from the soft 
and collinear gluon radiation which is independent of masses |5|,|§, 

tf04^~l + ^(- 2 + y + ■■••) • (4) 

Hence an accurate prediction for the rate to 0(ot z s ) can be obtained by calculating the K factor 
in the large M t limit and multiplying this by the lowest order rate computed using the full mass 
dependence. 

Along this line, we compute the QCD corrections to the processes e + e~ — > tth° and pp — > 
tih° in the M h /M t -> and high energy limits, which we expect to be reasonably accurate 



2 



from our experience with the gg — > h° process. In these limits the Higgs boson can be treated 
as a parton bremsstrahlung off a heavy quark. Section 2 computes the distribution of Higgs 
bosons in a heavy quark, ft-+hi to leading order in 77=, ^ and to 0(a s ) in the QCD 
coupling. In Section 3, we use the results of Section 2 to compute the QCD corrections to the 
process e + e~ — > tih°. We then combine the O(a^) corrections to pp — > it with the 0(a s ) 
corrections to ft->h t° compute the K factor for pp — > ith° in the high energy and large M t 
limits. Just as is the case for the gg h° process, this K factor can be combined with the 
lowest order result computed with the full mass dependences to obtain an accurate estimate of 
the radiatively corrected cross section. Section 4 contains some conclusions. 



2 Higgs Boson Distribution Function: The Effective Higgs 
Approximation 

2.1 Lowest Order 

For a Higgs boson which has a mass much lighter than all of the relevant energy scales in a 
given process, Mh « M t , the Higgs boson can be treated as a parton which radiates off of the 
top quark. This approach was first proposed by Ellis, Gaillard, and Nanopoulos ||, and we will 
refer to it as to the Effective Higgs Approximation (EH A). The process qq(gg) — > ith° can be 
thought of in this approximation as qq(gg) — > it, followed by the radiation of the Higgs boson 
from the final state heavy quark line. This procedure is guaranteed to correctly reproduce the 
collinear divergences associated with the emission of a massless Higgs boson. We expect that 
treating the Higgs boson as a massless parton should be a good approximation for the purpose 
of computing radiative corrections and the K factor. 

The distribution function, ft-,h, of a Higgs boson in a heavy quark, t, can be found from 
the diagram of Fig. ^| using techniques which were originally derived to find the distribution 
of photons in an electron |I0| . This approach has also been successfully used to compute the 



distribution of W bosons in a light quark [|TT|] . 

The spin- and color-averaged cross section for the scattering process te — > tX through 
t— channel Higgs (h°) exchange is, 



• I A(te -> tX) | 2 dT x , (5) 

where the kinematics are given in Fig. [|, vt and v e are the velocities of the incoming particles, 
and dTx is the Lorentz invariant phase space of the final state X. We factorize the amplitude 

as 
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Figure 2: The scattering process te — > tX can be used to extract the Higgs boson structure 
function ft-*h( x h), assuming that the amplitude is dominated by the t-channel Higgs pole. 



A(te -> tX) | 2 =| A{eh° -> X) | 2 ( 1 Y | A(t -> t/i°) | 2 . (6) 

V /b h 



This factorization is the crucial step in defining the Higgs distribution function and assumes 
that the dominant contribution comes from an on-shell Higgs boson. 
We define the energy carried by the Higgs boson to be 

E h = x h E , (7) 
and neglect all masses, except when they lead to a logarithmic singularity.^ Using 

a(eh°^X) = — i— dT x (27r) i 5 i (k+p e -P x )\A(eh°^X)\ 2 , (8) 

and d 3 p ' = (1 — Xh) 2 E 3 dxhdQ' we find 

E 2 f f dVt 1 

a(te -+tX) = — dx h x h {l -x h ) j — | A(t th°) | 2 a{eh° -> X) . (9) 

We now define the Higgs boson structure function by the relation, 

a(te = ^ f t ^ h {x h )a(eh° -> . (10) 

This is a Lorentz invariant definition which can be used beyond the leading order. We use the 
notation f^ h (xh) to denote the lowest order result, 



1 It is straightforward to include masses in the analysis. However, experience with the effective W approxi- 
mation has shown that the mass corrections are typically numerically unimportant. 
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where A(t — > th°) = —ig t 5iju(p')u(p). 

We now consider a Higgs boson which couples to heavy quarks with the Yukawa interaction, 

C = -g t tth° . (12) 

In the Standard Model, gt = — , with v = 246 GeV. Our calculation is, however, valid in any 
model where « M t , such clS db SUSY model. 

Performing the integral of Eq. (|TTD , we obtain our primary result, 



The dominant numerical contribution comes from the (1 — Xh)/xh term, i.e. from the infrared 
part of the Higgs distribution function. In fact, we can also reproduce the first term in Eq. ([Hp 
if we calculate, in the eikonal approximation, the bremsstrahlung of a soft Higgs from the final 
state of the e + e~ — > ti process. This provides us with an independent check of the leading 
behavior of the Higgs structure function we are going to use in the following. 

It is instructive to compute the neglected mass dependent terms using a different approach. 
The cross section for qq — > tth° has been computed many years ago [12| , 



ala(qq -> tth°) 
dx h 



27ns X 



, 4M, 2 - Mi 
x h + 2' • 



2 (AM? - Ml)(2M? - Ml) 8M, 



+ 



sx h 



log(^±|) 



4/3 



2M?\ /AM? — Ml 



)}• 



(14) 



where Xh = 2E h / y/s is the exact equivalent of the quantity defined in Eq. (|?p and varies in the 
range between 



xf" = ^ and x™ ax 



AM? Ml 
1 L + ^L 

s s 



(15) 



while f3 is defined to be 



P 



(x 2 h - (x^ in f)(x^ ax -x h ) 
(x™ ax - x h + AMl/s) 



1/2 



(16) 
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Figure 3: Virtual (a) and real (b,c) QCD corrections to be computed in order to obtain the 
0(a s ) corrections to the structure function ft^h(xh) m the EHA. 

We can define the distribution of Higgs bosons in a heavy quark by the relation, 

a(qq tth°) = 2 J dx h a(qq -> tt)f t ^ h (x h ) , (17) 

where the factor of 2 reflects the fact that the Higgs boson can radiate from either quark. Using 
Eq. ( |17|) we find an alternate definition of the Higgs distribution function, 



ft^hM 



dcr(qq — > tth) 



-2a (qq — > tt)\ dxn 
Taking the limit s » M? » M% of Eq. QT^D, we find 



16iU " "V Mf I ' V x h 

Except for the argument of the logarithm, this agrees with Eq. ([13]). It is easy to convince 
ourselves that this discrepancy can only introduce a difference of the same order of magnitude 
as the effects we are neglecting and is therefore irrelevant in our approximation. We have 
explicitely checked that the use of ft^h{xh) in the form of Eq. (fl3|) or of Eq. (19) is numerically 



irrelevant for an intermediate mass Higgs boson and parton sub-energies above around 1 TeV. 
2.2 QCD Corrections 

The utility of the results of the previous section is that it is straightforward to compute the 
0(a s ) corrections to ft^h- The discussion here parallels that of Ref. JTB| ]. 

2.2.1 Virtual Corrections 

The virtual corrections include contributions from the vertex renormalization of Fig. |3]a, from 
the top quark wavefunction renormalization, and from the top quark mass renormalization. 
The virtual correction from Fig. |||a is given by, 



6 



A"(t^th°)^C F 

47T 



2\ e 



4(1 + 2e)IW(e) - 6 



+ 



(20) 



where we define = M\jM^ and CV = |. ^4.° (t — denotes the tree level amplitude for 
t — > th°, and can be directly read from Eq. ([12]). We have explicitely separated the infrared 
(e/ij) and ultraviolet (euv) divergences and we have defined 



Tie) 



-~1E=- ■ 

e e 



:2D 



The wave function renormalization (Z 2 ) contributes to the virtual amplitude, 



Z 2 A°(t -> th G 



-A°(t 
-A°(t 



th c 



th c 



a. 



^ 47T/i 
47T/i 2 



2\ e 



4tt V M? 



a. 



4tt V M? 



C h 



l + 2e)T uv (e) + 2(l + e)T IR (e) 
2 



+ 



+ 4 



:22i 



The sum of (|2|) and (|20| ) has both an UV pole and an IR pole. The UV pole in fact is cancelled 
by the renormalization of the mass term in the Higgs-fermion coupling (g t = M t /v), while the 
IR pole will be cancelled by the soft gluon bremsstrahlung contribution (see Fig. |||b and ||]c). To 
complete the calculation of the virtual part, we add to ( |22|) and (^) the mass renormalization, 
i.e. 



rW(e)(3 + 4e) 



and we obtain, 



A 



virt 



A°(t 



th°)^C F 



47T 



2 r M 



3 -«-.« N -, . M 2 



— 6 r h + 0(r 



(23) 



(24) 



Therefore, the virtual processes contribute to the QCD corrected Higgs distribution function, 



In the limit — > we have 



2 r m f 47r ^ 



-6 



f*7. + 0(r 2 ) 



7T 



(25) 



(26) 
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2.2.2 Real Corrections 

The gluon bremsstrahlung diagrams are shown in Fig. ||]b and |||c. In the soft gluon limit, the 
amplitude is given by, 

A (t _> thPg)^ = g s T"A°(t - th°) ( tlh. _ ElEi.) . ( 2 7) 

The contribution to the Higgs distribution function is found by integrating over the gluon phase 
space, 



(PS) g = I dE 9 I sin 1 - 26 9 1 d9 1 I sm- 2 *e 2 dB 2 , (21 



-8 pl-2e 

JL i 1--'' « / 

'0 '(27r) Jo JO 
where is the gluon energy. The integral of Eq. ( pq) is cut off at an arbitrary energy scale, 
5, and the dependence on <5 must vanish when the hard gluon terms are included. Performing 
the integrals of Eq. ( p8f) , the soft gluon bremsstrahlung term is given by 

/mTW = t»W^ ||r fc r( CjJI ) + A + 0(r 2 )} , (29) 

In Eq. Q2T3D, A is the sum of all the remaining finite contributions, 



- r ,lo g (4) + -log^— j+.tog^^^ 

ciu 1U Ml-/3„(«) 



" (2 rh) io [1 - «(1 - u)r h \ ' (30) 



M? \ 1/2 , „ / M 2 \ V2 



and we have denoted by 

Pp=( 1 -^) and ^=^-^J (31) 
the velocities of the incoming (p) and of the outgoing (p') t-quark respectively, while 

is the /9-term velocity for the 4— vector v = up + (1 — which we use in the integration 
of the interference term between the two real emission diagrams. In the limit — > and 

M t «E,E', 



1 f /E + E'\ (E'\ 1 



(33) 



We can see that in the sum of (|25|) and ( p9|) the residual IR divergences cancel. 
The final QCD corrected result is 



ft~h{x h ) = f^ h (x h )\l + ^C F \-l + ^^\og(l-x h )}} + f hard + 0(r h ) 

I 7T L X h J J 

= fl h (* h ){l + ^fUM} ■ (34) 

The hard gluon terms cannot be reproduced in this approximation and they would require the 
full calculation of the bremsstrahlung of a hard gluon in the final state of e + e~ — > tih or 
pp — > tih . Note that to leading order in there is no dependence on the gluon energy cut-off 
S in the soft bremsstrahlung contribution of Eq. ( f29|) and so f hard can contain only finite terms. 
In addition, there can be no 1/xh singularities in f hard and so we expect this term to be small. 
The emission of a hard gluon should have a very distinguishable experimental signature and 
should not affect the study of e + e~ — > tih or pp — > tih°. Therefore we do not include the 
study of hard gluon emission in this context and drop the f hard part of the structure function 

Figure ^ shows the lowest order and radiatively corrected results for 2 typical energy scales. 
It is clear that the dominant contribution is from the (1 — Xh)/xh term with little energy 
dependence and that the order 0(a s ) contributions are small. 



3 Cross Sections 



3.1 e + e- -> tth° 

As a first application of the formalism previously derived, we want now to estimate the impact 
of QCD corrections on the cross section for e + e~ — > tih . This will require us to use the QCD 
corrected distribution function derived in Eq. ( [34]) together with the NLO expression for the 
cross section for e + e~ — > ti. 

Let us start testing the validity of the EHA approach in the absence of QCD corrections. 
For this purpose, we evaluate and compare the cross section for e + e~ — > tih in the exact theory 
(i.e. the Standard Model) and in the EHA, for an intermediate mass Higgs boson. 

The cross section for e + e~ — ► tih has been first calculated in |TJ]] (photon-exchange con- 
tribution only) and then completed in [15[ (both photon and Z-exchange contributions). To 



make the comparison with the formalism explained in Section 2 easier, we will write the cross 
section in the form 

u(e+er tih") = f " dx k Me+e ] - ^ , (35) 

J 'min CtXh 
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Figure 4: The Higgs distribution function ft-+h{ x h) i n units of g 2 j{A^ 2 . 

where Xh = ZEyJ y/s and the integration bounds have been defined in Eq. fll5|) . The differential 
distribution in Eq. (|35|) can be written as 



da(e + e~ 



tth c 



N, 



<T0 



(4tt) 2 
V 2 + A 2 



Q2Q i + 2Q e Q t V e V t , (y e 2 + Al)(V t 2 + A*y 



1 - Aff/s 



+ 



A 2 t Y J G l + V t 2 {G^ + G { 



i=2 



1 - Mils) 2 

QeQtVeVt 

1-Mf/s 



Gp 



(36) 



where a = 47ia 2 /3s, a e is the QED fine structure constant, and Qi, Vi and A4 (i = e, t) denote 
the electromagnetic and weak couplings of the electron and of the top quark respectively, 



V 



2P 3L -4Q t s 2 w 



9 P 

As — 2hL 



4s w c w ' As w c w 

with 7^ = ±1/2 being the weak isospin of the left-handed fermions and s\ ¥ = \ — = 0.23. 



(37) 
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The coefficients Gi and G2 describe the radiation of the Higgs boson off the top quark (both 
photon and Z boson exchange), while G3, . . . , G§ represent the radiation of the Higgs boson off 
the Z boson. Explicit expressions for the coefficients Gi, . . . , Gq are given in the Appendix. 
As already observed in Ref. [Tj|, the most relevant contributions are those in which the 



Higgs boson is emitted from a top quark leg, i.e. those proportional to G\ and G2 in Eq. (|36f) . 
This is further confirmed when we calculate the same process in the EHA and we see that it 
is numerically consistent with the exact result to within a factor of 1.5 — 1.8 for 50 < < 
150 GeV and y/s = 1 TeV (the agreement is clearly improved when we go to higher values of 
y/S)- 

In fact, according to Section 2, the cross section for e + e — > tth° is evaluated in the EHA as 
the convolution of the Higgs boson distribution function with the cross section for e + e~ — > tt, 
i.e. 

x max 

a(e + e- -> tth°) EHA = 2 / dx h f^ h (x h )a(e + e~ -> tt) , (38) 

/-.rain 
° X h 

and therefore only the emission of a Higgs boson from a top quark is taken into account. 

In the absence of QCD corrections, f^_, h {xh) is given by Eq. (|i~3"D, while cx(e + e~ — > tt) can be 

easily calculated and reads |l6|,|r7[ 



a(e + e -> tt) = ^ (l + °vv + Pt^AA , 



(39) 



where \i t = 2M t /y/s, /5 t 2 = (1 — /i 2 ), and ayv an d &AA correspond respectively to the product 
of two vector or two axial currents (the interference between the two gives zero upon angular 
integration) , 



N c a 



QlQ] + 2Q e Q t V e V t 



s-M% 



+ V t 2 (V e 2 + A 2 e ) 



Mi 



oaa = N c a A 2 t (V e 2 + A 2 e 



Mi 



(40) 



using the same notation we introduced before. We notice that cr(e + e~ — > tt) does not depend 
on Xh and the associated Higgs production, in the EHA, is obtained by simply multiplying 
cr(e + e~ — >• tt) by a prefactor, as can be seen from Eq. (RSI). 

In Figure |5| we compare the exact cross section of Eq. (|3^) to the EHA one of Eq. (p8|), in terms 
of their ratio 



R 



a (e 1 e 



a e T e 



-> tthP) 
tth°) 



EHA 



EXACT 



(41) 
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Figure 5: Ratio of the lowest order cross section computed in the EHA for e + e — > tth° to the 
same quantity computed exactly, for intermediate mass Higgs bosons at y/s = 1 TeV. 

As we can see, the EHA reproduces the exact cross section to a good level of approximation 
and we are entitled to use it in the estimate of the impact of QCD corrections on the associated 
production of an intermediate mass Higgs boson in ti events. 

Within the EHA formalism, we can parameterize the effect of QCD corrections as 



(7(6 e 



tth c 



dxy 



h fL h M f i + ^-fUM 



7T 



(42) 



where f^ h is the Higgs boson distribution function without QCD corrections, as in Eq. (|T3|) , 
while the exact expression of fl^^x^) can be derived from Eq. fl34|). Moreover, we have 
schematically represented the QCD corrected cross section a(e + e~ — > ti) in the form 



a[e + e 



it) = (T H o 



71 



(43) 



where <7 is given in Eq. (p9|). 
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Numerical expressions for a 1 at various energies can be found in the literature [17] and we refer 
to these numbers in our discussion. The impact of QCD corrections on the prediction of the 
cross section for e + e~ — > tih° is described by the K-factor, 



In the EHA, using Eq. 



58T) and Eq. 



a e 1 e 



*tth 



ith°) 



NLO 



a(e+e- -> tth°) LO ' 
, the K-factor can be written as 



(44) 



K(p) t 



*tth 



'i:f dx h ft h M (i + ^fuw)] 






-.max 

J x min dXfiff._^fa (X/x) 


6° 



1 + 



+ 



7T 



CP' 



x rnax 

J x min dXhfl—thiXh) 



(45) 



Using the results of Ref. JL7], we can estimate that a 1 / a Q at \/s = 1 TeV (taking fi = y/s 
and a s (s) = 0.088), while the second term in parenthesis is roughly —4 for < 180 GeV. 
Therefore K(s) ~ .94 and the impact of QCD corrections is estimated in the EHA to be very 
mild for e + e~ initial states. 



3.2 pp 



tth° 



The QCD corrections to the process pp — > tih° can be estimated using the EHA and the results 
of Ref. J7|,|| for the NLO corrections to pp — > it. The NLO cross section for it production at a 
hadron collider can be conveniently parameterized at the parton level as, 



a(ij -> itX) 



a 2 s (p) 



y -2 { fij(p) + 4 ™*G") [fjjip) + fl(p) log(^ 



where ij corresponds to the gg, qq, qg, and qg initial states and 



4M t 2 



(46) 



(47) 



with s the parton sub-energy. Analytic results for the /£(/>) and fy{p) and numerical param- 
eterizations of fh(p) can be found in Ref. [[/]]. 

In the EHA, the cross section for pp — > ith° production at lowest order is then, 



a°(pp -> tth°) 



dxidx 2 Fi(xi)Fj(x 2 )fii(p) / dx h ft_+ h (x h ), 



'.i 



(4S 
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where F^x) are the parton distribution functions and X\ and X2 are the momentum fractions 
carried by the incoming partons. Again, the over-all factor of 2 reflects the fact that the Higgs 
boson can be radiated off either top quark leg. The parton sub-energy is s = Xix 2 s. Note that 
because of the dependence on s in f t -,h(xh), the effect of the Higgs emission is not a simple 
prefactor as was the case for e + e~, (although the energy dependence is very small). 

At lowest order, the process pp — > ith° includes both the qq and gg initial states. As was 
the case for e + e~, the contribution from the qq initial state is well approximated at lowest 
order in the intermediate mass Higgs region by the EHA calculation. However, the process 
gg — > tth° also contains the t -channel emission of a Higgs boson which is not included in 
the EHA approximation and so the EHA is a much poorer approximation to the lowest order 
hadronic cross section than it is to the e + e~ cross section.^ We will therefore use the EHA only 
to calculate the K factor. The best estimate of the rate for pp — > tih° will then be obtained 
by using the exact calculation for the lowest order rate and then multiplying by the K factor 
obtained using the EHA. 

The cross section to NLO can easily be found using the EHA, 



cr\pp -> tth C 



Ml 



dx 1 dx 2 ^2F i (x 1 )F j (x 2 ) / dx h fl^h( x h){fij(p) 



71 



fUM + (^ 2 ) /i(p) + 4(p)io, 



EL 



(49) 



The K factor is then given as usual by, 



^{pp -> tth° 



*tth 



(50) 



181 It 



a°(pp -> tth°) 

In Fig. ^|, we show the K factor using structure functions derived by Morfin and Tung, 
is clear from this figure that using the lowest order cross section with the 2— loop evolution of 
a s (n) overestimates the size of the corrections. (This is also true for the pp — > h° process. ||]) 
The K factor varies between around 1.2 and 1.5 for the intermediate mass Higgs boson. Since 
the dominant production mechanism is gluon fusion, the exact value of the K factor is sensitive 
to the choice of structure functions. We have found that changing the set of structure functions 
induces a 10 — 20 % uncertainty in the K factor shown in Fig. M. 



2 The agreement between the exact calculation and the EHA calculation of pp — > tth° at the LHC can be 
improved by using Eq. ( |l9| ) instead of Eq. ( ]l3| ) for f^ h - Even so, the EHA overestimates the exact cross 
section for Mh ^150 GeV by about a factor of 2. The K factors obtained using Eq. ( |l3| ) or Eq ([l9|), along with 
Eq. (|34|), are, however, almost identical. 
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Figure 6: K factor for pp — > tth° production at the LHC, y/s = 14 TeV, for /i = M t , using 
Morfin-Tung [18] structure functions. The solid curve uses the 1— loop evolution of and 
lowest order structure functions for cr°, and the 2— loop evolution of a s (/i) and NLO structure 
functions for a 1 . The dashed curve uses the 2— loop evolution of a s (fi) and the LO structure 
functions for both a and a 1 . 

4 Conclusions 

We have computed the next-to-leading order QCD corrections to e + e~ — > tih° and to pp — > tih° 
in the high energy limit and the limit in which Mh << M t . We expect this to be a reasonable 
approximation in the intermediate mass region. The corrections to e + e~ — > tih° are small and 
can safely be neglected. The corrections to pp — > tih°, however, increase the rate at the LHC 
by a factor of between 1.2 to 1.5 for < 180 GeV, although the numbers are sensitive to the 
choice of structure functions. Within the context of the EHA, the bulk of these corrections can 
be identified with the 0(a^) corrections to the pp — > ti sub-process. The significant size of the 
corrections underscores the need for a complete calculation. 
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A Coefficients for a(e + e —> tth°) 



As we have seen in Eq. fl36|), the differential cross section for e + e~ — > tih° can be expressed 
in terms of six coefficients G%, . . . , Gq. Two of them, G\ and G2 account for the emission of a 
Higgs from a top quark and are explicitely given by 



G x = 

s 2 [p 2 -x 2 h ) x h 



, 2g2t k ( ~ 4/3 (4M 2 - M 2 h ) (2M 2 + s) x h + (51) 
[p 2 - x 2 ) a;?, V 

p 2 - xf) (l6M t 4 + 2M^ - 2M^sx/, + s 2 x 2 h - AM 2 (3M^ - 2s - 2sx h )) log 



^ +J3 



G 2 = . 2g * r (^a* (-96M 4 + 24M 2 M 2 - s (-M 2 + s + so*) (-(3 2 + x 2 h )) + 

s 2 (p 2 -xl)x h \ V • . V J J 



2 (p 2 - x 2 h ) (24M 4 + 2 {Ml - M 2 h sx h ) + M 2 {-UM 2 h + 12sx h + sx 2 h )) log 



x h +_$ 
x h - P / 

while the other four coefficients, G3, . . . , G$ describe the emission of a Higgs from the Z-boson 
and can be written in the following form, 

Gs = M2( "'fff - 2 m + 12M* z + 2M 2 z sx 2 h + s 2 (-l + x h )x 2 h - 

M| (M^ - M z + s- sx h ) 

M 2 (8M 2 +s(-A + 4x h + x 2 ))) , 

G 4 = \ 9 } 2 Ml - 2 U8Ml + 12M 2 h -s(-2A + p 2 + 2Ax h -3x 2 h )) , (52) 

6 (Ml — MJ + s — sx h ) v v ' ' 

G = = Mz , ( - Ml+ ^M\ + s ( - 1 + Xk)) ^ + '» - ' + »>)) + 

2 (M 2 (M 2 - 3M| + 8 - sx h ) + M 2 (-4M 2 + 12M| + sx 2 )) log 



x h + P 



x h - P t 

8g t g z M t M z ( * . 2 2 , fx h + p 

G « = 3 (-A g + M I + 3 (-l + ^) ^ S + (4M ' " " 4 + 2S " ^ '° g 1^, 

In both Eq. flBTD and Eq. (|52|) , (ft and denote the Yukawa couplings of the top quark 
(g t — M-tjv) and of the Z boson (gz = Mz/v) respectively, while the quantity j3 is defined in 
Eq. fll6[) . Ref. [pjl presents the result for da / dxxdx2, where x\ and x 2 are the fractional energies 



of the top quarks. After making a change of variables and integrating over | x\ — x 2 |, their 
expressions for the Gi agree with ours. 
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